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We simulate binary soft disk systems in two dimensions, and investigate how the dynamics slow as 
the area fraction is increased toward the glass transition. The "fragility" quantifies how sensitively 
the relaxation time scale depends on the area fraction, and the fragility strongly depends on the 
composition of the mixture. We confirm prior results for mixtures of particles with similar sizes, 
where the ability to form small crystalline regions correlates with fragility. However, for mixtures 
with particle size ratios above 1.4, we find that the fragility is not correlated with structural ordering, 
but rather with the spatial distribution of large particles. The large particles have slower motion than 
the small particles, and act as confining "walls" which slow the motion of nearby small particles. The 
rearrangement of these confining structures governs the lifetime of dynamical heterogeneity, that is, 
how long local regions exhibit anomalously fast or slow behavior. The strength of the confinement 
effect is correlated with the fragility and also influences the aging behavior of glassy systems. 



PACS numbers: 61.20.Ja, 64.60.My, 64.70.pv, 81.05. Kf 



I. INTRODUCTION 



Many liquids can form glasses if they are cooled 
rapidly, and glassy materials have technological applica- 
tions such as optical fibers and plastics The origin 
of the glass transition is still unclear, despite the scien- 
tific and technological interests. Much work has exam- 
ined the dynamical properties of materials near the glass 
transition. Those studies revealed several important fea- 
tures of supercooled liquids and glasses. For example, 
upon approaching the glass transition, the structural re- 
laxation time (tq.) increases by several orders of magni- 
tude without a corresponding growing static correlation 
length 0, The rate of this increase of Tq, is called 
"fragility" and depends on the material [H-IH . For fragile 
glass-formers, Tq, steeply increases for a small decrease 
in temperature, while for "strong" glass-formers, the in- 
crease of Ta requires a larger decrease in temperature. 
The relaxation time is related to the viscosity, and thus 
the fragility is an important factor in the ease of pro- 
cessing glass-forming materials. Typically it is desirable 
to mold a glass-forming material with the viscosity held 
within a certain range; for fragile materials, this may cor- 
respond to a rcstrictively narrow range of temperature. 

Another common observation of materials close to the 
glass transition is that they often have a broad distri- 
bution of local mobility: some regions in a sample relax 
much faster than other regions. Eventually, molecules 
in those regions exchange their dynamics, that is, fast 
regions become slow and vice versa. This is termed "dy- 
namical heterogeneity" [sl-IT^ . It is known that the char- 
acteristic lifetime of these dynamically heterogeneous re- 
gions TjjH is longer than Tq and the strong divergence 
of tdh near a glass transition temperature suggests that 
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heterogeneous dynamics are relevant for understanding 
the glass transition Q. 

A third common feature of materials in the glassy state 
is that their properties evolve with time, such as the dif- 
fusivity of molecules or dielectric susceptibility. This phe- 
nomenon is termed "aging" [T5l420j . Despite the chang- 
ing properties, no clear structural changes have been 
seen; this can be true even if, for example, diffusivity 
slows by several orders of magnitude [l^ [2^ . 

While all of these properties have been known for some 
time and carefully characterized by experiments, the ori- 
gins of many of them are unclear. To understand the 
origins, numerical simulations with simple intermolecu- 
lar interactions are used to study factors controlling the 
dynamics. These model systems are useful as clear under- 
standing is hindered by the complexity of real materials 
[2l| - [25j . For the fragility, some simulations show that liq- 
uids become less fragile when the polydispersity increases 
(or a larger size ratio for binary mixtures is used) (22| - [23 | . 
This has been explained as due to ordering of the sample, 
which becomes frustrated in polydisperse samples. For 
example, in two-dimensional (2D) systems, small regions 
with hexagonal order can form which correspond with 
slower dynamics (larger values of Tq,), and thus increasing 
the polydispersity frustrates formation of these ordered 
regions and diminishes the fragility [l^. Furthermore, 
it is reported that particle mobility in those ordered re- 
gions are slower than that in randomly structured regions 
and it suggests that dynamical heterogeneity is also in- 
fluenced by local structure [l^ [25l - l33| . 

Those simulations used the polydispersity as a small 
perturbation frustrating the ordering, that is, with either 
a small polydispersity or a binary system with the parti- 
cle size ratio close to 1. However, dynamics in those sit- 
uations are quite different from highly polydisperse sam- 
ples. Binary mixtures have two control parameters: the 
size ratio and the volume fraction ratio of the two com- 
ponents. These lead to complex phase diagrams [H, HH 
and potentially emergent dynamical properties such as 
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an effective depletion interaction between the large par- 
ticles in binary hard sphere systems (36l.[37|. Our interest 
is in dense amorphous phases at intermediate size ratios. 

In this article, we simulate the dynamics of binary soft 
disk mixture systems with large size ratios and find that 
the glass transition in large size ratio binary systems can 
be quite different from that of systems with smaller size 
differences. We study the fragility of binary soft disk 
mixtures with various size ratios and area fraction ratios. 
Local ordering is less significant for large size ratio sys- 
tems. Instead, we see a molecular crowding effect from 
the large particles [H, Our data suggest that the 
large particles act as confining walls for the smaller par- 
ticles, and that confinement effects increase the fragility 
of such systems by suppressing dynamical heterogeneity. 
We also investigate aging in large size ratio systems, and 
again seen an influence of confinement effects due to the 
large particles. Overall, our results suggest that confine- 
ment effects are crucial to describing the dynamics of the 
glass transition in binary mixtures with large size ratios 
between the two components. 

II. METHOD 




FIG. 1: (Color online.) The orientation pair correlation as a 
function of distance. The straight line corresponds to 
the behavior expected for hexatic phases. A sample with 
size ratio aL/crs = l-2 and area fraction ratio 4'l/4's ~ 0.5 
at = 0.66 does not have hexatic order since g^i^R) / g{R) de- 
cays faster than the power of —1/4 (the lower curve, marked 
"Glass"). In contrast, the upper curve (marked "Hexatic") 
shows the correlation for (t_l/(ts = 1.15, 4>l/4>s ~ 0-5, and = 
0.66, which has long range orientational ordering. 



We perform two-dimensional Brownian dynamics sim- 
ulations for binary mixtures composed of large (V) and 
small (5) soft particles. These simulations are meant 
to mimic the colloidal glass transition. For the colloidal 
glass transition, the key control parameter is the vol- 
ume fraction 0, [Til, S^] , and so in our simulations the 
chief control parameter is the area fraction The re- 
sults of Brownian Dynamics simulations are similar to 
those of molecular dynamics simulations in dense systems 
[H [23. liTI - liat . The particles interact via the purely re- 
pulsive Weeks-Chandler- Andersen potential |4J|; C/y = 



4e[(cr,y/r)i2 _ (cr.y/r)6 -f 1/4] for r < 2^^^a,j, otherwise 
Uij = 0, where (Ty = {ai + aj)/2 and i,j G {L, S}. We fix 
ksT/e = 0.04, so the total area fraction (j> is our control 
parameter to approach the glass transition. The mass ra- 
tio is mL/fns={o'L/ (^sY ^'^'^ the length is normalized by 
as- The total number of particles is iV = Nl+Ns = 1024 
(or 4096) where Nl and Ns are the number of large and 
small species, respectively. We generate initial conditions 
by simulating at 4> = 0.10 for a long time and expanding 
the disk sizes to increase (j>. We confirm that our results 
do not show initial condition dependence and our results 
are also independent on time below (j) = 0.66, that is, 
there is no aging observed. Thus, we consider our binary 
systems are well mixed and ergodic for (f> < 0.66. 

For certain area fractions for a 2D monodisperse sam- 
ple, the hexatic order phase can be found, which has ori- 
entational order but no long-range positional order [4^ . 
For binary samples with large size ratios, there is ap- 
parently little ordering of either type, while for size ra- 
tios close to 1, hexatic phases are still subtle to verify. 
Hence, we calculate an orientation pair correlation func- 
tion ge(i?)/g(i?) where i? is a normalized distance from a 



particle. ge{R) / g{R) decays with the power of R^^''^ in a 
hexatic phase for a monodisperse sample [45| . Thus, we 
classify our samples as liquid (or glassy) when the power 
of the decay is faster than R~^^'^ [Fig.[l]. We find no hex- 
atic phases for any area fractions for systems with size 
ratios cri/(7s=1.2, 1.25, 1.5, 1.75, 2, 2.5 and 3 and area 
fraction ratios 4>l/ 4'S=NL<^\/Ns(j'g=Q-?>j 0.75, 1, 1.5 and 
2; these states are indicated by the circles in Fig. [2ja). 
We confirm that the decay rates of ge{R)/g{R) in those 
simulation conditions arc quicker than i?^/"* and thus we 
arc studying liquids (or glasses at higher (f>). Wc note 
that it is known that several crystalline structures ex- 
ist in three-dimensional (3D) binary sphere suspensions 
[3^ . l46j . In none of our samples do we observe large 
hexagonally ordered patches for our simulations, and the 
small patches that sometimes appear are transient. This 
is discussed further below. 

When the size ratio is large, a depletion force should be 
generated [11] . This has been studied before in dilute sus- 
pensions of large and small spheres, and manifests itself 
as an effective attractive force between the large spheres 
[13] ■ We cannot rule out the existence of the depletion 
force in our simulations, which may affect the structure 
of the large particles, although at large area fractions the 
depletion force may be less relevant. Indeed, wc observe 
pairs of adjacent large particles which neighbors for long 
periods of time, but these neighbors separate eventually 
(see Movie SI, EPAPS materials to be submitted), and in 
general all of the dynamics are slow in our glassy exper- 
iments. It's unclear if depletion introduces still slower 
dynamics, or if this is merely part of the overall slow 
behavior. At least, we note that our systems are quite 
different from gels with strong attractive force. 
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FIG. 2: (color online). The contour plot of (a) fragility index 
D, (b) the growth rate of hexagonal order d{%pQ)/d<j), (c) the 
growth rate of dynamical correlation length 9^/90 and (d) the 
dynamical correlation around a large particle O in a plane of 
{(Jl/o-s, 0l/<^s)- All those figures except (a) are obtained at 
(j>=0.66. The shading in each figure is darker when each value 
is small, thus, the darker regions correspond to (a) fragile 
liquids, (b) low growth of hexagonal order, (c) low growth of 
dynamical correlation length and (d) minimal correlations of 
motion between a large particle and its neighbors. The circles 
in (a) are the simulated points and the specific states A, B, 
C, D, E and F are labeled. 



III. RESULTS 

A. Fragility: behavior as (f> increases 

We obtain the relaxation time Tq, from the self-part 
of the intermediate scattering function for all particles, 
which is given by F{k,t)^j^ X]j(expjfc • [rj{t) - ^^(0)])' 



where fj is the position vector of particle j, () indicates 
a time average and k is the wave vector, is deter- 
mined when F{kp,Ta)=l/e where kp corresponds to the 
wave number of the first peak of the structure factor. 
The (j> dependence of Ta is well fitted by Vogel-Fulcher 
function substituting for 1/T: Tq^tq exp[D(j>/ ((/)o — (/))], 
where D is the fragility index and 4>o is the area fraction 
of the ideal glass transition [see Fig. ^a.)] 0. Fragile 
liquids have smaller values of D. For example, D ^ 4 
for triphenyl phosphite which is one of the most frag- 
ile liquids, while for the less fragile liquid butyronitrile, 
30 [ig. For our simulations, 0.4 < Z? < 1.0, smaller 
than the molecular liquids. The difference may be due 
to using the density as the control parameter rather than 
temperature. For comparison, we examined the data of 
Refs. [H, \^ which used light scattering to study the 
colloidal glass transition as a function of volume frac- 
tion. From their data, we find D = 0.497 ± 0.002. An- 
other example is that the fragility index of glycerol in 
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FIG. 3: (Color online.) The area fraction (f> dependence of 
(a) the relaxation time Ta, (b) average hexagonal ordering 
("i/je) and (c) the dynamic correlation length ^ at states A, B, 
C, D, E and F indicated in Fig. [^Ja). The solid lines in (a) 
are fitting lines with the Vogel-Fulcher function at each state. 
Filled and open symbols correspond to fragile liquids and less 
fragile liquids, respectively. 



an isothermal experiment is smaller than that measured 
in an isobaric experiment [Hoj . We expect our observed 
qualitative dependence of the fragility on the system pa- 
rameters to still be revealing. 

It is worth noting that for our data, we also compute 
the value of D by calculating F{k,t) for only small par- 
ticles (or large particles) and we obtain almost the same 
values for D and (j>o. We do not see separate glass tran- 
sitions for the two particle species, which has been seen 
in simulations with large size ratios and equal particle 
numbers (thus volume fraction ratios 4'l/4>s much larger 
than ours) [H, HI]. In such cases, the large particles do 
not seem to interact as directly with the small particles 
but rather have their own glass transition, and then the 
small particles move in the interstices between the large 
particles. In our simulations, the volume fraction ratios 
are somewhat close to 1 , and so the large particles always 
"see" the small particles and the alpha relaxation time 
scales for both particle species have similar (p dependence. 

Figure ^a) shows the contour plot of the fragility D 
in a {phjasi 4'l/4>s) plane. The circles indicate the 
states simulated. This figure shows that the fragility 
is a non-monotonic function of the size ratio or the 
area fraction ratio. For example, for 4'l/4's = 2.0, 
the fragility has values D = 1.05,0.82,0.45 for states 
A, C, and E, but then increases slightly to 0.49 for the 
state to the right of E. Likewise, for 4'L/<t>S = 0.5, the 
fragility behaves non-monotonically with increasing size 
ratio: D = 0.64, 0.48, 0.77 for states B, D, and F. Com- 
paring states A and B, or C and D, suggests that increas- 
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ing the number of large particles increases the fragility 
index Z?, but states E and F disprove this trend. 

We find that fragile liquids (small D) have small 0o, 
the area fraction where Tq appears to diverge. The re- 
lationship between fragility and the divergence point of 
Tq is observed at glass-forming liquids [sil [52j and our 
results are consistent with them. On the other hand, the 
existence of 0o for molecular liquids is still discussed and 
not well established [H, H^l • We are not sure what de- 
termines 4>o and why fragility is related to the divergence 
point. Below, we focus on microscopic properties such as 
particle mobility and local arrangement. 

Prior work observed that two-dimensional systems can 
form small hexagonally ordered regions, and the mobil- 
ity of particles is diminished within these regions. More 
fragile liquids are observed to have large growth rates of 
the size of these regions with respect to <t> [22 2j]. To 
check this we study the </> dependence of hexagonal order 
for our samples. We use the local hexatic order param- 
eter described as i/'e = ;r" I Y^m=i e*^^'" | where Uj is 
the number of nearest neighbors for particle j and O-j-^ is 
the angle of the relative vector r,„ — rj with respect to 
the X axis. "06=1 means that a hexagonal arrangement 
is formed around particle j, while V'e^O corresponds to a 
non-hexagonal arrangement. We then consider (ipe), the 
time and particle average of '0g. We compare the 4> de- 
pendence of {iPq) with that of Ta [Fig.|31^a) and (b)]. Both 
particle sizes are similar at states A and B, olI^^s ~ 1-2, 
and here both Tq and ('06 } increase dramatically as in- 
creases, suggesting they are indeed related as seen in prior 
work. However, for states C and E, (-06) stays nearly 
constant with increasing 0, while Tq, grows rapidly. The 
system slows without significant hexagonal ordering. 

Next, we compare the fragility with the growth rate of 
hexagonal ordering 9(V'6)/30 since fragility corresponds 
to d log Ta/dcj). (We calculate all derivatives of a quantity 
X with respect to as dX/dcj) = [X(0)-X(0- A0)]/A0 
where we choose A0 = 0.01. In the results below, we 
compare the behavior of the samples at a fixed = 0.66. 
However, our results show similar trends when compare 
samples at fixed r^.) Figure [21[b) shows the contour 
plot of d{tp(i)/d(j) at = 0.66 in a plane of {<tl/<Js, 
4'l/4>s)- The growth of hexagonal order is small at large 
ctl/cts and (j)L/4>s (upper right region). This behavior is 
expected since hexagonal ordering should be frustrated 
with increasing Ohj <^s- Figure|31Ja) shows Z? as a function 
of 9(^/'6)/90. We observe two distinct behaviors. For sim- 
ilar particle size systems {il/us < 1.4, triangles), more 
fragile liquids (smaller D) have a larger dependence of 
hexagonal order on 0, in agreement with prior work (2^ . 
In contrast, large size ratio systems {(Tl/cts > 1.4, cir- 
cles) show less correlation between the growth of hexag- 
onal order and the fragility index. 

Cooperative motion of groups of particles is a com- 
mon phenomenon as the glass transition is approached 
[lOl - fTsl . [55I [56j , and this behavior is thought to be more 
common in fragile glasses. Figure[51Ja) and (c) show snap- 
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FIG. 4: (Color online.) The fragility index D as a function 
of (a) d{Tp6)/dct> and (b) 9^/90 at 0=0.66. The triangles 
correspond to systems where size ratio is close to 1 (a_L/as < 
1.4) and the circles correspond to large size ratios (> 1.4). 
The lines are added as guides to the eye. 

shots of the systems at = 0.66 at state B and state E, 
where particles are colored based on Ar^, and groups 
of highly mobile particles are seen (darker colors). To 
define displacements, here (and for the subsequent anal- 
ysis in this work) we focus on the time scale At* for 
which the no n- Gaussian parameter is the largest. The 
non-Gaussian parameter is defined as 

-'^*>-^-' (" 

with displacements Ar measured over the lag time At, 
and the factor of 3/5 chosen so that 0:2 = for a Gaussian 
distribution f57| . a2 is larger when the tails of the dis- 
tribution become broad. Prior work identified the time 
scale At* that maximizes a2 as related to cage rearrange- 
ments [1, 0, [III . In our current work, we computed a2 
separately for the large and small particles, finding sim- 
ilar values for At*. For our analysis, we will use At* 
based on a2 calculated for the small particles, and our 
results are not sensitive to this choice. 

Examining our data, we observe clusters with cooper- 
ative motion in our systems, some of which arc circled in 
Fig. [5l[a,c). To look for the connection between the co- 
operative motion and fragility in our sample, we need to 
characterize the cooperative motion. We compute a cor- 
relation function described as S'(i?) = (Ari-Arj)/(|Arp), 
where R is the distance between particles i and .7 , and 
Afi is the displacement of particle i at At=At* [55[ . We 
find that S{R) shows exponential decay with R, which 
was previously observed in experiments [56j and simula- 
tions This exponential decay yields a decay length 
^, which we plot in Fig. as a function of for our 
six representative states. While all samples have similar 
short ranged cooperative motion at « 0.60, we see a 
variety of behavior as the glass transition is approached. 

Figure El^c) shows little relationship with the magni- 
tude of ^ and the fragility, so we focus on the growth 
rate of ^, d£,/dcf>. This is sensible since the fragility in- 
dex D relates to the growth of Ta as increases. Fig- 
ure [21^c) shows the contour plot of d^/dcf) in a {(Jl/cS: 
4>l/4's) plane. This plot has a rough qualitative similar- 
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FIG. 5: (Color online.) (a) and (b) are snapshots of the sys- 
tem at state B, </> = 0.66; (c) and (d) are at state E, (/> = 0.66. 
The simulation box is four times as large as those snapshots in 
each side. In (a) and (c), particles are colored based on mobil- 
ity Ar|, with darker colors indicating more mobile particles. 
The darkest color corresponds to Ar^=0.25a^ . The displace- 
ment time scale At is chosen to maximize the non-Gaussian 
parameter, and is At = 5 x 10* for (a) and 2 x 10* for (c). 
In (b) and (d), particles are colored based on the hexagonal 
order parameter -i/jg, with the darkest color corresponding to 
"05=0.8. The outlined regions are guides to the eye. 



ity to the fragility [Fig. [2ja)]. Figure HJ^b) more directly 
shows that D is related to d£,/dcf>, with distinct behav- 
iors for (Jl/o's < 1.4 systems and aL/<ys > 1-4 systems. 
In Ohl'^s < 1.4 (triangles), fragile liquids have larger in- 
crease of ^ with respect to and it is natural that the 
increase of Tq is related to this. This is also suggeste d by 
experimental studies of colloidal suspensions [hL l56l | . 
In contrast, the opposite relationship is seen for the large 
size ratio states [circles in Fig. lU^b)]. The most fragile 
states (small D) correspond with those where ^ grows 
least as the glass transition is approached. This unusual 
behavior seems to be a key result for understanding the 
fragility of large size ratio binary mixtures. 

To further understand the relation between ^ and 
fragility in large size ratio systems, we consider the dy- 
namical difference between the two particle species. We 
compute the mean square displacement (Ar^) of large 
and small species separately at a variety of states, and 
find that the large particles are always significantly slower 
than the small particles. We next consider how the mo- 
tion of small and large particles is coupled. To quantify 
this, we compute the correlation between the directions 
of motion of a large particle and its neighboring small 
particles as = {-^ cos 9ij) where j indicates a large 
particle, i are the nearest neighbor particles for particle 
j, rii is the number of these neighbors, 9ij is the an- 



gle between Ari and Afj, Afj is the displacement of j 
particle at At=At*, and the angle brackets indicate a 
time average and an average over all large particles j. 
= 1 indicates that particles around a large particle 
move in the same direction as the large particle on this 
time scale, while = means that their movements are 
uncorrclated with the large particle. Figure HJd) shows 
the contour plot of 8 in the plane of (ctl/cts, (I)lI4>s)- 
ranges from 0.7 (upper left corner) to 0.4 (lower right 
corner). Cooperative motion decreases as the size ratio 
increases, and as the number of large particles decreases. 

If large particles move slower and small particles move 
independently of the large ones [the lighter-shaded re- 
gion in Fig. md)], this suggests that the large particles 
act as slow-moving walls within the large size ratio sys- 
tems [1^ . The small particles are trapped in pores 
between the large particles [see Fig. [5jc)]. In confined 
geometries, can dramatically change compared to an 
unconfincd system at the same temperature and density 
[58l - l60| . Confined systems with free surfaces typically 
have smaller values of T n, w hile those with rigid walls 
have larger values of Tq, [61|. Our data suggest the lat- 
ter situation is relevant for our large size ratio systems, 
that crowding due to the large particles slows the motion 
of the small particles. This has been seen in prior work 
[38ll39| and here we examine this behavior in more detail. 

Here, we suggest the fragility of large size ratio sys- 
tems may be connected to the "strength" of confinement 
effects. Less fragile liquids may have more mobile walls, 
such as state C with a large value of Q (implying small 
and large particles move together). Or, less fragile liquids 
may have a larger spacing between the large particles, 
such as state F, with a relatively small value of (j)L/4's] 
here, they are less confined. In contrast, the more fragile 
states D and E have small values of 9 and smaller dis- 
tances between large particles. These systems are ultra- 
confined, where the spacing between large particles is of 
order ^ or even smaller, thus limiting 9^/90. 



B. Influences of particle mobility at constant (j) 

We also investigate the relationship between local 
structure and local mobility in our binary systems. Ac- 
cording to prior work, the mobility of particles decreases 
when the particles are in hexagonally ordered regions (for 
2D simulations) [22], [25l. which is also hinted at in 3D 
colloidal experiments [2^. Thus, we focus on the spatial 
distribution of mobility and hexagonal structure. Figure 
[SJa) shows a snapshot of the system in state B at (/)=0.66 
with the darker colors indicating particles with larger val- 
ues of Ar|, and Fig.[5l^b) shows the same snapshot color- 
ing the particles by their value of ?/;g . The circled regions 
in Fig. [Sl[a,b) show that mobile regions correspond to 
regions with less ordering. We calculate the Pearson cor- 
relation coefficient C(Ar|,?/;g), finding C = —0.13, sup- 
porting the idea that mobility is slightly anticorrelatcd 
with hexagonal ordering, consistent with prior work. 
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FIG. 6: (Color online.) (a) Schematic showing a group of 
small particles temporarily confined within a region bounded 
by large particles. Each small particle is assigned a value 
Rw which is its closest distance to any large particle (the 
large particles act as "walls"). Within each region, the local 
maximum of identifies the "confinement pore size" Rc 
for that region. Thus, Rc is the maximum radius of a circle 
that fits within the bounded region, (b) A contour plot for 
the mean square displacement (Ar^) (at fixed At = At") 
as a function of the effective pore size Rc where the given 
particle is located, and the distance of that particle from the 
nearest "wall" R^, for state E with (p = 0.66. Darker shading 
corresponds with high mobility. The motion of particles is 
faster in larger pores (large Tic) and when further from the 
walls (larger R^ for a constant -Rc). The length scales are in 
terms of the small particle diameter. For these data, i?* ~ 3.6; 
only 10% of regions have Rc > Rc- 



For large size ratio systems, we see little hexagonal 
structure on average [curves E and F in Fig. [31[b)], but 
this does not preclude the possibility that locally there 
may be hexagonal ordering which influences the dynam- 
ics. To check this, we compare the local mobility with 
local structure. Figure [5jc) and (d) show Ar^ and -0g 
for state E, with a much larger size ratio, and here there 
is no correspondence (C = 0.01). Overall, for the large 
size ratio samples (ctl/cts > 1-4), we never observe any 
ordered structures at any area fraction. However, we can- 
not rule out the possibility that there is subtle ordering 
that might be present and influencing the dynamics. 

Next, we consider the relationship between the confine- 
ment effects and mobility in large size ratio systems. The 
confinement effects are composed of both confinement 
size effects and confinement surface effects. For confined 
colloidal suspensions, particle motion was slower in more 



confined spaces, but there was not a strong infiucnce 
from the confining surfaces [59l |. Numerical simulations 
show that the particles move slowly near rough walls, 
and quickly near smooth walls [s^ . We thus wish to dis- 
tinguish between finite size effects and interfacial effects. 
First, we define clusters of large particles as those large 
particles separated by a distance less than IAcl, and 
these clusters form "walls" surrounding small particles. 
In some cases, a connected cluster of large particles com- 
pletely surrounds a group of small particles, as sketched 
in Fig. |6l Within such a region, we compute the distance 
Rw from each small particle to the nearest "wall" parti- 
cle. The maximum value of Rw within the confined region 
defines Rc, the effective "confinement pore size." R^ and 
Rc are indicated in Fig. |6l[a) . Ru, is calculated per small 
particle, and Rc per pore; both of these are functions of 
time. The dependence of specific particles' behaviors on 
Rw should give insight into interfacial effects, and the de- 
pendence of pore-averaged behavior on Rc should give a 
separate insight into finite size effects, although of course 
these two effects are likely both present simultaneously. 

Figure [6lb) shows a contour plot of (r^) as a function 
of a (i?t„, Rc) plane at state E and cf) ~ 0.66. The re- 
sults arc located only at the lower right of the graph as 
Rc > Rw fi'om our definitions. The darker region at the 
upper right corresponds to high (r^), showing that the 
mobility increases inside large pores (large Rc) and far 
from walls (large Rw)- Given the incommensurate sizes 
of the large and small particles, it is difficulty for small 
particles to pack well near the "walls" [see Fig. [5jc)] and 
so not surprisingly our results arc consistent with simu- 
lated rough walls [s^ . There is also a slight gradient of 
increasing mobility as a function of pore size Rc for fixed 
Rw < 1.5, indicating that there is a finite size effect in 
addition to an interfacial effect. That is, smaller pores 
(smaller Rc) have more particles close to the pore walls, 
and thus experience stronger interfacial effects, but the 
data indicate that the influence of the interface on adja- 
cent particles is less within large pores. 

To obtain further evidence for the relationship between 
the conflnement effects and dynamical heterogeneity, we 
also investigate the temporal relationship between local 
structures and dynamical heterogeneity. We calculate 
the intermediate function F{k, At) for small and large 
particles separately at fixed k = 2Tr/ as for small parti- 
cles (thick solid (black) line in Fig.[71[a)) and k = 27r/(Ti 
for large particles (thick dashed (blue) line in Fig. [71[a)). 
The structural relaxation time scales for small and large 
particles, ts and tl, are set by F{k,T) = 1/e. We next 
consider all of the confined regions and the distribution 
of region sizes Rc- We determine the distribution of all 
pore sizes Rc (taken over the entire simulation run), and 
find the threshold size i?* for the top 10% of this distri- 
bution. For each pore at each time, we define Wc{t) = 1 
if Rc{t) > R* and Wc{t) = otherwise. Typical values 
of Rc range from 2.8 to 5.6. The temporal correlations 
of the regions are given by {WdAt) ■ Wc{0)} / {Wc{0)^} , 
plotted as the thin solid (red) line in Fig. ^a). The 
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FIG. 7: (Color online.) (a) The time dependence of F{k,At) 
for small particles (thick solid (black) line; k — 2h/(ts) and 
for large particles (thick dashed (blue) line; k = 2n/(TL) at 
state E and (j) = 0.66. The thin solid (red) line indicates the 
breakup of the largest confined regions, and the thin dashed 
(green) line indicates when the most mobile particles become 
less mobile, (b) The typical structural relaxation time rg for 
small particles (filled circles), the structural relaxation time 
tl for large particles (open circles), the confinement lifetime 
Tc (diamonds), the dynamical heterogeneity lifetime tdh (tri- 
angles), and the peak time for the non-Gaussian parameter 
At* (squares). The data are for state E. (c) Five time scales 
{tl, Tc, Tdh, At* , and Thex) normalized by ts as a function 
of ctl/cs at 4'l/4's = 2.0 at (/) = 0.66. r^ex is the lifetime of 
hexagonal order and corresponds to the filled hexagon sym- 
bols. (These data go through state points A, C, and E; see 
Fig. W) (d) Five time scales {tl, Tc, tdh. At* , and Thex) 
normalized by rg as a function of (jl/cts at 4>l/4'S ~ 0.5 at 
(j) = 0.66. (These data go through state points B, D, and F.) 
There are only a few confinement regions when ahjos < 1.5, 
preventing a clear determination of Tc, so those curves are 
clipped in (c) and (d). 



typical lifetime of large regions is given as Tc where the 
correlation drops to 1/e. Similarly, for each particle we 
define a parameter Wdh which is equal to 1 if the parti- 
cle's displacement Ar at that time is within the top 10% 
of the displacement distribution, and otherwise. The 
correlation (WoHiAt) ■ M^d//(0))/(M^dh(0)2) is plotted 
as the thin dashed (green) line in Fig. ^b), and toh is 
defined by the 1/e time again; this is the time scale over 
which particles exchange between being fast and slow, as 
mentioned in the Introduction. 

Figure [7Jb) shows the 4> dependence of these four time 
scales [ts, tl,Tc, tdh and At*). The fastest time scale is 
Ts (filled circles) , and the large particles are much slower 
(tl, open circles). What is more notable is that tl ^ Tc, 
in other words, large particle motions relate to the relax- 
ation of confined regions. This is further evidence that 
the large particles form walls. Furthermore, Tc ~ tdh is 
observed, connecting the confinement-induced dynamics 
with the dynamically fast particles. Faster particles ex- 
change identities with slower particles when the confining 
walls rearrange. 

If the spatial dynamical heterogeneity is actually in- 
duced by the confinement effect, we would expect to 
see Tl K, Tc ~ tdh, and the correspondence should be 
strongest for large size ratio systems. Figure[71Jc) and (d) 
show these time scales normalized by T5 as a function of 
the size ratio ul/cfs at 4>l/4'S = 2.0 and 0.5, respectively; 
the data are for </) = 0.66. Indeed, for aj^jas > 1-5, we 
find Tl is similar to Tc and t^h- Not surprisingly, the 
relaxation time scale Tc for the confinement effect is gov- 
erned by the relaxation time tl for the large particles 
which define the pores; more significantly, the lifetime of 
dynamically heterogeneous regions is also connected to 
the behavior of the large particles. 

On the other hand, we find /ts 1 when tJi /as ~^ 
1, which is to be expected. In these cases, tdh remains 
large, showing that here the dynamical heterogeneity is 
less influenced by the relaxation of the large particles. 
However, r^j^f is same order as the lifetime of hexagonal 
order, as seen in Fig. [ZIc,d) by comparing the triangles 
(tdh) with the filled hexagons (r/ie^). This new time 
scale, Thex, is defined in a similar way as the other time 
scales: particles with 'ipe > 0.8 are considered hexagonal, 
and the correlation time for having ipQ > 0.8 is Thex- This 
time scale is only relevant for samples with reasonable 
amoimts of hexagonal order [compare Fig.[5l^b,d)], and so 
is only shown for samples with small size ratios. (This is 
why it is not shown in Fig.[7|b) which has (Tl/ct^ = 2.5.) 
It is further evidence that local hexagonal order influ- 
ences the dynamics in similar size ratio systems. In these 
cases, the confinement effect is more likely due to hexago- 
nal regions composed of both particle species, rather than 
networks formed by only the large particles, and thus Tc 
becomes ill-defined and we do not show it in Fig. [7ljc,d). 

The competition between hexagonal ordering and 
crowding due to large particles likely accounts for the 
cases in Fig.[7ljc,d) where tdh > tl- For example, states 
with ai^jas = 1.5 still have regions of hexagonal or- 
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dering composed of both particle sizes. In these cases 
Q is smah [see Fig. [l{d)] and the confinement effect is 
weak. Another example where tdh > is the state 
with aL/as = 3.0, 0l/</'s = 0.5 [Fig. EJd)]. The large 
particles are scarce, but result in a strong confinement 
influence; however, the more numerous small particles 
can themselves form hexagonal patches, influencing their 
mobility. We expect that for 3D glass-formers, local crys- 
talline order is much less significant, and so confinement 
effects would more strongly determine tjjh in all cases. 

Furthermore, we consider the relationship among At* 
and other time scales (Fig. [7l[c) and (d)). As a reminder. 
At* is the time scale at which the small particle displace- 
ment distribution is the most non-Gaussian. This time 
scale is used to define the mobility of particles and so is 
part of the definition of a "mobile" particle and thus tdh- 
At* > Ts shows that small particles move appreciable 
distances during the time At*. That is, the non-Gaussian 
displacements are over significant distances, enough to re- 
lax the small particle structure, and this becomes more 
true at larger size ratios, although this motion is all local- 
ized within a region defined by nearby large particles. In 
all cases tdh > At* showing that slow and fast regions do 
not change identities with each particle rearrangement, 
but rather take longer times to change. Particles may re- 
arrange several times confined within a large pore (sev- 
eral At*) before the pore rearranges and the particles 
change their mobility. An important caveat is that these 
observations are for the average behavior of particles, and 
we are not implying that the connections between pore 
sizes and dynamics are strongly deterministic. Nonethe- 
less the hierarchy of time scales in Fig. [7] is suggestive of 
nontrivial connections between the spatial arrangements 
of the large particles and the long-lived lifetimes of dy- 
namically unusual regions, both fast and slow. 



C. Aging 

Next we investigate aging dynamics in binary samples 
with large size ratios. We study state E at = 0.72 where 
it is in a glassy state {(j)g ^ 0.68 for this state). Figure 
IHlJa) shows (Ar^) for small and large species separately at 
waiting times t.^, = 10^, 10"', 10^ and 10^. At short time 
scales (At < 5000), (Ar^) is independent of t^. Particles 
move within cages formed by their nearest neighbor par- 
ticles. At longer time scales At > 5000, (Ar^) increases 
as the cages rearrange and allow particles to move to new 
positions. This upturn in (Ar^) decreases with increasing 
tu,, indicating the aging of the system, similar to what 
has been seen in experiments |18l - [20| . Figure [DJb) shows 
the tu, dependence of (Ar^) at fixed At = 10^. Though 
the results at small t^ {tw < At) are hard to interpret as 
the dynamics change during At, we can observe the clear 
temporal change of (Ar^) for both particle species. For 
"old" systems, the plateau of (Ar^) extends over a large 
range of time scales, with the plateau height correspond- 
ing to the cage size. Within our uncertainty, the large and 
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FIG. 8: (a) (Ar^) as a function of lag time At for small 
particles (solid lines) and for large particles (dashed lines) at 
t„ = 10^, 10", 10^ and 10® at state E and </> = 0.72. The 
mobility of both particle species decreases with respect to 
aging time, (b) The dark black lines correspond to (Ar^) at 
At — 10^ as a function of aging time tm for small particles 
(solid) and large particles (dashed). The smooth curves are 
fits with a stretched exponential, and we find that the decay 
time of large particles is slightly shorter than that of small 
particles. The gray line show the waiting time dependence 
of Q, the correlation between the displacements of a large 
particle and its surrounding small particles. The choice of 
At = 10'' is arbitrary although panel (a) shows that this is a 
reasonable choice to capture the slowing dynamics. 



small particles age almost at the same rates, as suggested 
by the similar shapes of (Ar^) curves (Fig. [5Ja)). Again, 
these observations are consistent with experiments in bi- 
nary colloidal glasses j20| . 

We wish to know a reason for the slowing dynamics 
with respect to the waiting time. First, we confirm that 
the overall structure is unchanged with age: we compute 
the pair correlation function g{R) at t.^ ~ 1000 and 10^ 
and cannot observe any difference between them, similar 
to prior observations in simulations [6^ and experiments 
[iMIl^j- g{.R) is a spatial average over the whole system, 
so we next consider the relation between local structure 
and local mobility. We use the confinement pore size Rc 
to characterize the local structure of the particles within 
a given confining region. Figure [9] shows (Ar^) and the 
probability distribution P(i?c), both as functions of Rc at 
t^ ^ 1000 and 10*^. We find that P{Rc) does not change 
at all and this is further evidence that the structure does 
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FIG. 9: (a) The probability distribution P{Rc) of confine- 
ment pore sizes at tn, — 10^ (black line) and t-w = 10^ (super- 
imposed white line), showing the distributions are essentially 
identical. This is for state E with (j) = 0.72. (b) Mean squared 
displacements (Ar^) for fixed At = 10^ as a function of the 
confinement pore size Rc at t^, ~ 10^ (circles) and tw = 10*' 
(squares), for the same state as (a). The upward trend (faster 
motion for larger pores) is similar to what is seen in equili- 
brated liquids [Fig. EJb)] . For these data At = 10^ to capture 
time scales over which the dynamics age [see Fig. |8|b)]. 



not change. However, the mobility decreases with age 
with the dependence on Rc relatively unchanged other 
than the amplitude. This implies that the mobility de- 
crease is not due to local structural changes, but rather 
an average slowing of the whole system [l9[ . 

Next, we focus on the importance of confinement which 
strongly influences the dynamics of our equilibrated large 
size ratio binary systems. As noted above, the confine- 
ment strength depends on the confinement size Rc and 
the effective rigidity of the walls, that is, 0. In aging 
systems, P{Rc) does not change with respect to tt„, but 
Q could depend on t^i since 8 is a dynamical property 
rather than a structural property. We compute 8(iu,) 
at fixed A< = 10 ^, shown as the gray line in Fig. [HJb). 
The behavior of Q looks similar to the mobility change 
of both particles. Figure [TU] shows a scatter plot of the 
mean square displacement (Ar^) as a function of 6, for 
all large particles and all waiting times t^, . We can clearly 
see the correlation between the mobility and Q. When 
@ decreases, the cooperative motion between small and 
large particles is less, in other words, the rigidness of walls 
increases. This result implies that confinement effects be- 
come stronger during aging and it may help explain the 
slowing down of the mobility. However, we don't know 
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FIG. 10: (Color online.) Scatter plot of the mean square 
displacement (Ar'^) as a function of at state E and (j) = 
0.72, using At = 10^. Each point corresponds to a different 
tm, the time since the start of the simulation. The solid line is 
a least squares fit. The data indicate that mobility is linked 
to the correlation of the motion between large particles and 
their nearest neighbors. 

why 8 decreases as the sample ages. 

IV. CONCLUSION 

We have examined the glass transition in binary mix- 
tures with a large size ratio, finding results that are dis- 
tinct from binary mixtures with smaller size ratios. Sys- 
tems with smaller size ratios are often studied, and the 
utility of using two particle sizes in those cases is to frus- 
trate the packing and prevent crystallization. Crystals 
are also frustrated in our simulations with large size ra- 
tios, and in addition we find several new results. First, 
we have investigated fragility of binary systems. The 
fragility of binary systems with size ratios close to 1 
((Tl/(T5 < 1.4) is related to the growth of ordered regions: 
more fragile liquids show dramatic increases in hexago- 
nal order as the glass transition is approached |22h24| . 
However, systems with larger size ratios do not show this 
relation. Both types of systems have fragilities which 
are related to the growth of a dynamical length scale ^, 
although the sign of this correlation is opposite for the 
small and large size ratio systems. 

The data show that in large size ratio systems, large 
particles act as quasi- immobile walls which confine the 
small particles and slow the dynamics overall. The large 
particles define regions with a range of sizes. Small parti- 
cles in large regions find it easier to move, even if within 
that region they may be adjacent to a large particle at 
the boundary of the region. Small particles in smaller 
regions are much less mobile. In addition to these fi- 
nite size effects, there are interfacial effects: small parti- 
cles near large particles move slower than those farther 
away. Those results are what is often seen in experi- 
ments and simulations of confined supercooled liquids. 
Furthermore, as is often observed in simulations and ex- 
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pcrimcnts, we find some particles are unusually mobile. 
Our new observation is that the length of time which 
these particles stay unusually mobile is connected to the 
lifetime of the confinement effect and thus the relaxation 
time of the large particles. 

Finally, we also investigate aging dynamics in large size 
ratio systems. Below the glass transition, the mobility 
decreases with respect to the waiting time, though we 
cannot observe any structural change. We find that in 
"younger" systems that the motion of large particles is 
correlated with the motion of their neighbors, but that 
in "older" systems this correlation is markedly smaller. 
This correlation (or lack of it) relates to the confinement 
effect, suggesting that the large particles become more 
rigid confiners in older samples. 

It is important to note that simulations of softer par- 
ticles with a charged (Yukawa) potential find results dif- 
ferent from ours, pointing out that our results are not 
completely generalizable. In simulations with a size ratio 
1 : 5, the large particles crystallized [1^. In these cases, 
the large particles did not rearrange but rather moved 
on their lattice sites, and small particles could only move 
by diffusive hopping motions between the crystal inter- 
stices. This is in contrast to our simulations where the 
large particles are always able to rearrange (albeit more 
slowly than the small particles). 

Overall, our results suggest that in binary soft sphere 



systems, the effect of the large particles to induce finite 
size effects within the sample play an interesting role in 
the dynamics. Two relevant variables are the finite sizes 
of regions between large particles, and the effective rigid- 
ity of the large particles. While our simulation studies 
2D systems, the results are similar to prior observations 
in 3D binary colloidal experiments [20|, |62| with moder- 
ate size ratios. While the effects are easiest to see with 
large size ratios, one implication of our simulations is that 
these effects may be relevant although less obvious in bi- 
nary systems of smaller size ratios. Indeed, one of the 
goals of our simulations was to understand the effect of 
structure by using systems where structural heterogene- 
ity is more obvious. These results may also have implica- 
tions for studies of nanocomposites, where inclusions into 
polymer glasses can dramatically affect the properties of 
materials [g^, [fi^ . 
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